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Abstract: We investigate the theoretical description of the central exclusive production 
process, hi + hi — > h\ + X + hi. Taking Higgs production as an example, we sum loga- 
rithmically enhanced corrections appearing in the perturbation series to all orders in the 
strong coupling. Our results agree with those originally presented by Khoze, Martin and 
Ryskin except that the scale appearing in the Sudakov factor, /x = 0.62\/I, should be re- 
placed with /i = y/I, where \f§ is the invariant mass of the centrally produced system. We 
confirm this result using a fixed-order calculation and show that the replacement leads to 
approximately a factor 2 suppression in the cross-section for central system masses in the 
range 100-500 GeV. 
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1. Introduction 

At hadron colliders, in events producing high transverse momentum particles in the central 
rapidity region, the colliding particles usually break up. However, in a small fraction of 
events the colliding hadrons remain intact and scatter through small angles. This type of 
production is known as central exclusive production (CEP): 

h(pi) + h 2 (p 2 ) -> hi (pi) © X © h 2 {p' 2 ) , (1.1) 

where the © denote rapidity gaps between the outgoing hadrons and the central system X. 

The outgoing hadron momenta may be parametrised in terms of the momentum frac- 
tions each transfers to X, Xi, and their transverse momenta, p', , : 

p'f = (1 - xOpjf + -^_p^ + pf ± (1.2) 
j# = (1 - x a K + tt^vK + (1.3) 
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with s denoting the centre-of-mass energy squared. The CEP kinematics are defined as 

^ < Xi < 1 . (1.4) 

X{S 

If the outgoing hadron momenta are measured, by adding detectors far down the beam- 
pipe, it is possible to reconstruct the four- momentum of the central system X. In the CEP 
kinematic regime, ( |l.4|) , the central system's rapidity, y, and invariant mass squared, s, are 
given approximately by 

s « X1X2S , (1.5) 



2 V^2 

In addition, the process possesses a P-even, C-even, J z = selection rule [1] (the origin 
of which we shall discuss subsequently). Thus CEP offers a method to measure both the 
mass of X [2] (with a resolution of ~ 2 GeV per event [3]) and its spin-parity properties [4]. 
Photon pairs [5], di-jets [6] and \c particles [7] produced via the CEP mechanism have now 
been observed at the Tevatron and there are groups within both the ATLAS and CMS 
collaborations actively seeking to observe these events at the LHC [3]. 

For a massive enough central system, the calculation of this type of process may be 
performed in perturbative QCD. At lowest order in the strong coupling, two gluons in 
a colour singlet state are exchanged between the hadrons. Crucial however to a correct 
theoretical description of the process is the inclusion of logarithmically enhanced terms, 
appearing at all orders in the QCD perturbation series. 

In this paper we investigate the theoretical description, in perturbative QCD, of the 
CEP process, taking Higgs production as an example. We begin in section ^ by introducing 
central exclusive production and describing the calculation by the Durham group. In 
section ^ we discuss central exclusive Higgs production at lowest order, highlighting some 
of the general features of the amplitude. Then, in section |], we sum large logarithms 
appearing at all orders in the strong coupling. Following this analysis, we see that the 
Durham result requires modification and in sections || and |6| we describe two fixed-order 
calculations that confirm this result. In section 0, we assess the phenomenological impact 
of our results. Finally, in section ||, we conclude. 

2. The Durham model 

The calculation of this process by the Durham group is represented schematically in figure |l|. 
The protons exchange a two gluon system, which must be in a colour singlet state in order 
that the protons remain intact. Two of the gluons then fuse to produce the central system, 
X. The cross-section is assumed to factorise in the following way [8,9]: 

dsdy dpf ± dpf ± ^ dady^^~^^' 

The transverse momenta of the final-state protons, p'^, are assumed to be distributed 
according to a Gaussian, with the slope parameter, 6 = 4 GeV~ 2 [9]. The factor S 2 , 
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Figure 1: Schematic form of the CEP amplitude. 



known as the soft survival factor, accounts for the suppression of the cross-section due 
to the requirement that soft interactions between the incoming protons do not spoil the 
exclusivety of the process [10]. In general S 2 is expected to depend on the kinematics of 
the process, however it is common practice to set it to a constant value, corresponding 
to the average over the forward detector acceptance of the final-state proton transverse 
momenta. In any case, the factorization of S 2 is not expected within QCD, although there 
are good reasons to suppose that it holds to a reasonable approximation. We shall not 
consider the matter of soft survival any further in this paper and defer instead to the 
recent literature [11—16]. 

The partonic cross-section, a, is related to the matrix element for two on-shell gluons 
to produce the central system as 

da(gg -+ X) = ^ \M{gg -> X) | 2 dPS x (2.2) 

As 

where dPSx is the phase-space of the final state, X, and 

M(gg -+X) = £ £ Ka 2 S XlX2 M a x \%(gg -»• X) (2.3) 

a\d2 A1A2 

with M.'x^ the amplitude for two on-shell gluons, with colours a, and helicities Aj, to fuse 
to produce X. Note that the averages are carried out at the amplitude level, in contrast 
to an inclusive partonic cross-section. In addition, the fact that the gluons have equal 
helicities (Ai = A2) is the origin of the J z = selection rule. 
Lastly, the effective luminosity is given by 

dC 1 ( vr fdQ] 

ds dy 

The f g are skewed, unintegrated, gluon distribution functions. Due to the kinematics of 
the process the amplitude is dominated by the region x\ <C x\ and in this regime these 
distributions may be related to the conventional, integrated, gluon density [8,17]: 

f g (x, x', Ql,n 2 ) « R 9 q^q2 (VT(Q±,^9(x, Ql)) • (2.5) 



^2^/^/ 9 (^^i,Qi,M 2 )/ fl (^2,4,Qi,M 2 )) • (2-4) 
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The factor R„ is given by 



Rn 



zt ( x x. r\l 
n 9\2 I 2 ; *** J 

xg(x,Q\) 



(2.6) 



and accounts for the skewed effect (H g is the skewed gluon distribution, see for exam- 
ple [18]). R g is approximately equal to 1.2(1.4) at the LHC(Tevatron) 1 [9, 19]. The f g 
distributions also include a Sudakov factor [17,20]: 

4 dfcj a s (kp r 1 -* 

IQ± k\ 2vr J 

where 

k± 



T(Q±,n) = exp f - J 



dz 



A 



k± + (jl 
H = 0.62Vs . 



(2.7) 



(2.8) 
(2.9) 



The Sudakov factor resums logarithmically enhanced soft and collinear virtual corrections 
and accounts for the fact that real radiation from the process is forbidden. 

The claim is that this expression resums logarithms in s/Q\, to next-to-leading log- 
arithmic accuracy. That is, it takes into account all terms of order a™ ln m (s/Qj_), with 
m = 2n,2n — 1. This requires a precise specification of both the lower limit on the k± 
integral in equation (2/7) and the cutoff on the z integral as z — > 1. Note that the upper 
cutoff on the kj_ integral corresponds to non-collinear hard radiation and as such there is 
no logarithm associated with this region. Thus, to next-to-leading logarithmic accuracy, 
only its order of magnitude is required. 

The lower cutoff on the k\ integral must be of the order of Q^, since radiation of 
a much lower transverse momentum would not be able to resolve the exchanged colour 
singlet system, the size of which is of order l/|Qj_|. To extract the precise value, the 
Durham group use the fact that this region, with k± ~ \Q±\ and the momentum fraction 
integral producing a logarithm, may be described within the BFKL framework [21-24] (see 
for example [25]). The BFKL summation of the momentum fraction logarithms implies 
the following replacement [26]: 



d 2 k ± 



d 2 fc ± 



1 



Q 2 



k 2 ± + (Q± - k ± f 



d 2 k ± 



Q 



(2.10) 



Thus determining the lower limit. We shall discuss this point in more detail in section 4.1 



Having specified the lower limit in this way, the Durham group fix the cutoff on the 
z integral, which they state is due to wide angle soft gluon radiation [27], by considering 
the cross-section for two on-shell gluons to fuse to produce a Higgs plus one additional real 
gluon, which they then argue, thanks to unitarity, may be used to imply the form of the 
virtual corrections making up the CEP Sudakov factor. To be more specific, they consider 



^or a LHC running at 14 TeV. 
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the region Q\ <^ <^ m 2 H , leaving the k± integral of the radiated gluon unevaluated but 
integrating numerically over its polar angle. The result is then fit to a function of the form 
a + 61n(m#/(2fc_L)), where a and b are /^-independent constants. This approach gives the 
following result for the cross-section [27] 

The factor of 11/12 here is the usual component of the f3- function coming from the z — > 1 
finite pieces of the P gg {z) splitting kernel. It is stated that the factor of 0.212 is due to the 
region of wide angle soft gluon emission. If this is the case, it may be absorbed into the 
logarithmic term (which is also due to soft emission). Equivalently, it may be written in 
terms of a momentum fraction integral 

a(gg -+ Hg) oc / / zP gg (z) dz (2.12) 



where A is given by equations ( \^g ) and ( |2.9[ ), Unitarity then guarantees that the soft 
cut-off in the real and virtual corrections must be identical, allowing one to infer the form 



of the z-integral cut-off in the Sudakov factor, (2.7). 



Our finding is that this result is not correct. Specifically, rather than equations (| 



and (2.9) we find instead A = k^/mu- In the next sections we will describe our evidence 
for this assertion: firstly in the form of an all-orders approximation to the CEP amplitude; 
secondly in an explicit next-to-leading order calculation of the relevant virtual corrections 
and thirdly, via a calculation of the gg — > gH matrix element. 



3. Lowest order Higgs production 

We begin our investigation of the Durham model by computing (in Feynman gauge) the 
lowest order amplitude for two quarks of different flavour to scatter into two quarks and a 
Higgs, Alo- We take the quarks as scattering in the forward direction: 

^ = (l-xiK, (3.1) 
p'? = (l-x 2 )p$. (3.2) 

The dependence of the amplitude on the outgoing hadron transverse momenta in equa- 



tion (|2.1|) is a non-perturbative effect and so we do not expect to be sensitive to it here. 
We will also make frequent use of the high energy limit: keeping only terms not suppressed 
by an inverse power of the centre-of-mass energy. In this limit, there are four lowest order 
diagrams which contribute, as shown in figure ^. In addition, we work in the effective 
theory in which the top quark has been integrated out [28-30]. This generates a point-like 
coupling of the Higgs to gluons (see appendix |A|) . 

For small Xi the amplitude is dominated by the region in which the exchanged gluons 
are soft and we may therefore take the gluons as coupling to the quark lines via eikonal 
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Pi, (71 



Pi, (71 
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Pi, (71 



Pi, 172 




p 2 ,<7 2 




Pi, (72 



Pi, (71 



P2, (72 



(b) 



P-,,(T 2 




(c) (d) 
Figure 2: The lowest order diagrams contributing to q + q' —> q © H © q', in the high energy limit. 



vertices. Thus, for example, the contribution to the amplitude of graph 2(a) , with the 
colour singlet contribution projected out, takes the form 



2^-1 4. 



N 2 



(2tt) 4 D 



(3.3) 



where 



£> _1 = [Q 2 + *e] [(Q ~ Pif + »e] \{Q + P2? + »e] 
x [(Q - x lPl ) 2 + ie][(Q + x 2 p2? + ie] 

and the Higgs vertex factor, H^ v , has the form (see appendix |A|) 

H^(h,k 2 ;n) = -iC^ih-h gT ~ KK) ■ 



(3.4) 



(3.5) 



Working in the centre-of-mass frame of p\ and p 2 , with p\ defining the z-axis, we may 
evaluate the Q + and Q~ integrals for each graph. To leading power in the high-energy 
limit we obtain 



M2A = . jPtol = = 



A LO = A 



2(b) 



dQl 



(3.6) 
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where A is a cutoff imposed to define the integral and 



Aq(jj) = -S aia >S na >C^{ji)g 



4 S TpCF 



(3.7) 



7T N 

Some comments are in order at this point. The first is that the amplitude is infrared 
power-divergent as A — > 0; this is to be expected, since we are dealing with a fixed-order 
expansion of a parton-level matrix element. As we shall see, this divergence will ultimately 



be tamed by a Sudakov factor (see also the footnote at the beginning of section 4.1). The 
second comment is that the contribution to the S-matrix (iA) is pure imaginary. This is 
expected for colour singlet exchange in the high energy limit, based on arguments from 
Regge theory [25]. Thirdly, note that it is perhaps not surprising that diagram 2(b)| gives a 
different result from all of the other diagrams, since it is the only one in which Q is pinched 
in the Glauber /Coulomb region [31,32]. 

Since its imaginary part dominates, we could equally well have used the Cutkosky 
rules [33] to evaluate the amplitude. In fact, it is this approach we shall use to evaluate the 
next-to-leading order corrections in section ||. The structure of the lowest order calculation 
will give us some hint as to how we may simplify the next-to-leading order corrections and 
so we detail it here. 

Recall that the Cutkosky rules implement the unitarity relation, which states that if 
we sum over the cuts of a graph, Gc, this sum is proportional to the imaginary part of the 
graph without the cut, G: 



J2°c = 23m(-»G) . 



(3.8) 



c 



For a cut to be allowed it must satisfy two conditions: the invariant-mass of the four- 
momenta on each side of the cut and the total energy flowing from left to right across the 




(a) 



(b) 



(c) 





(d) (e) 
Figure 3: The allowed cuts of the CEP amplitude. 
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Figure 4: The possible cuts of the lowest order diagrams. 



cut must both be positive. Given these conditions, the possible cuts of the general CEP 
amplitude are shown in figure [3| and the possible cuts of the lowest order diagrams are 
shown in figure III. 



However, it turns out that we need only consider the cuts [4(a)| and |4(e)| ; the remaining 
cuts all cancel amongst one another. The reason for this is that these diagrams represent 
a sum over soft gluon insertions onto the upper or lower quark lines. The fact that the soft 
gluon in all these cases is cut ensures that it is not in the Coulomb/Glauber region and we 
may therefore apply the soft approximation [32,34-36]. The result is that we may rewrite 
each sum of diagrams contributing to a given cut as a single diagram with the soft gluon 
connected to an eikonal line. This is shown explicitly in figure [| for the diagrams 4(b) 
and 4(g) (see for example [37] for the eikonal Feynman rules). The two diagrams differ 
only in the attachments of the gluons to the upper quark lines, which give the following 
expressions 



t b t a 



AM 



X 



oc 



(pi - Q) 
t a t b 

(p[ + Q) 2 



A* a 
2-PlPl 



U a 

PiVi 



+ 



~2 P1 -Q 

t a t b 
2pi-Q 



M a 

PxPi 



U a 

PiPi 



(3.9) 
(3.10) 



Here we took both gluons as soft, dropping all terms suppressed by Q or x\, which al- 
lowed us to make the eikonal approximation in the numerator and to set p[ ~ p± in 
the denominator of the second diagram. In addition, the on-shell delta- function gives 
|Q 2 | = \2x\pi ■ Q\ <C \2pi ■ Q\, allowing us to approximate the denominators as shown. The 
sum of diagrams is proportional to the commutator [t a ,i 6 ] or, put another way, in the 
eikonal diagram of figure | the upper quark line is connected to the lower part of the dia- 



-8- 



gram by a single gluon. The key point is that we are dealing with colour singlet exchange 
and such diagrams contribute only to the octet exchange part of the amplitude. We may 
therefore neglect them. 




Figure 5: Rewriting a sum of cuts in terms of an insertion onto an eikonal line. 



This situation generalises to all orders and so we need never consider the cuts 3(b) - |3(d) 



since they do not make a leading contribution to the colour singlet exchange amplitude. 
4. All orders Higgs production 

The inclusion of large logarithms appearing at all orders in the perturbation series is crucial 
if we are to obtain a good approximation to the CEP amplitude. Indeed, as we saw at 
lowest order, the result is divergent without them 2 . The logarithmically enhanced terms 
correspond to emissions that are collinear with external particles and/or soft. In the next 
sub-section we shall discuss the subset of these corrections corresponding to hard collinear 



emissions, deferring a treatment of soft emissions to section |4.2 . 

4.1 Hard collinear emissions 

4.1.1 Corrections that factorise into the pdfs 

We begin with the ansatz that the central exclusive Higgs production amplitude, -Acep, 

factorises at some scale, hf, much lower than all other scales in the problem. The amplitude 

2 Actually, the amplitude must be finite even without the inclusion of Sudakov effects, since as Q± — > 
the wavelength of the exchanged gluons becomes too large to resolve the colourless protons. Such an 
effect requires the introduction of hadronic wavefunctions and it has been studied in [38], using a simple 
non-perturbative model [39,40]. 
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Figure 6: Factorisation of the central exclusive Higgs production amplitude. 



may then be written (for small x\ and X2) as 



ao,ao 



where the kinematics are shown in figure |6| and we define 



P M = Pi + Pi , 



(4.1) 



(4.2) 
(4.3) 



The large logarithmic corrections to the coefficient function, C a ° a ° , may be computed 
by cutting off the parton transverse momenta at fj,p, which is equivalent to dimensional 
regularisation in the MS scheme [41]. The hard collinear corrections may be organised into 
two transverse momentum ordered chains as follows: 



c a °"°(£ ,T7,/4;£o,^/4) = E 



n=0 



E 

n=0 



E 
E 







dl n a s (l n ) 



4 l l ^ 



Ql dEa.(E 



fl' Tl fl' Tl — 1 



(0) 



E 



/? 47T 



&nQ"n — 1 



E 



(~ta n a n 

2 s-channcl 



(0) 



d6Kf f 



(4.4) 
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where Q± and Q± are the transverse momenta of the two s-channel ("rung") emissions 
closest to the hard scatter. The convolution symbol (g) indicates that these momenta are 
integrated over (e.g. see equation ( 4,16| )). The KlL are the analogue of the DGLAP 
splitting functions in the case of skewed kinematics. We have omitted their arguments for 
clarity but it is to be understood that 



K. 



(0) 



k (0) 



i-X 



They are given by [18] 

rqq 
'(0) 



Kfl(xi,x2\yi,y2) = C F 



- , &i 1 (xi,xi — yx) H — — fin(x 2 ,X2-y2) 



K q ^ ) (x 1 ,x 2 \yi,y2) = -y 



x\ - y\ " ' x 2 - y2 

+ i9i U (xi, -x 2 ,xi - y{] 

$112 O^i, -%2,xi - yx) + 2— — — , &xix(xi,-X2,xi - yx) 

2/12/2 



Kf^(xx,x 2 \yx,y2) = 2C F [(yx -2/2)^111(^1,-^2,^1 - yx) 



K 9 ^(xi,x 2 \yx,y2) = C A 



xx 

yi 



+XxX 2 ${xx(zii -X 2 , Xx ~ yx) 



'' 1 -^(xxjXx-yx) 



xx - yx 



+ X_2 
2/2 



-V 11 {X2,X2 - ?/2) 



+ 2 



X2 ~ 2/2 

xxx 2 xxyx + x 2 y2 q0 



xxX2 + yxm a0 , * 

+ 2 ^111(^1, -X2,xx - yx) 

2/12/2 

1 A) 



2/12/2 (xx + x 2 ) 2 
where the generalised step function is defined as 

ak 1 \ f 

"a\—a4 \ x li ■ ■ ■ 1 x j ) 



^ n ( Xl ,- X2 ) + [-^L+2)5(xx-yx) 



27ri ULx( x i K - 1 + i£ ) 



(4.5) 



(4.6) 



and the plus-distribution is given by 



-•d\\(xx,xx-yx) 



xx - yx 



' 1 -til^xxjxx-yx) 



xx - yx 



5(xx-yx) I dx'.-r^-VUx^x'x-yx) • 



(4.7) 



Note that our definitions of the kernels differ slightly from those of [18], in particular: 

1 



K {0) ~ K (0)\ [18] 



K (0) ~ 4^ K (0)\ [18] 



The collinear logarithms may now be absorbed into the pdfs: 

-4cep ^2sY^J^J diH a {^ m Ql)H- a {lmQl) C 2 i_ channel (£, r?; 77) 



(4.8) 



(4.9) 
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where 



H tt (£, v, Ql) = H a (S, 77; /4) + E / d ^ ^»(^o, V, /40 



n-l 



E -I1E 

n=l \ j=l a 3 J ^F "j 



'■Ui dlja s {lj) 



Ij 4vr 



-3 "(0) 



l (0) 



and likewise for iT a (£, 77; Qj_). 



?7 + Cn-l ??-Cn-i' 



(4.10) 



yip 



(a) 




Q 
(b) 




Q 

(c) (d) 
Figure 7: Diagrams contributing to Cf 9 s _ channcl . 



To study the coefficient function, we take {a, a} = {q, q} as an example. The diagrams 
contributing to Cf's-channel are shown in figure [j]. The contribution from figure 7(a) , for 
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example, may be written 



( ig 2 T F ^ f d 4 Q Tr [7-7^7" 



'ig 2 T F ^ f d 4 Q Tr [ 7 V$7o 



y m Z^J (2ir)*[Q* + ie][{y 1 p-Q) 2 +ie)[(y 2 p + Q) 2 + ie]J 

x \i{yw - Q) 2 (2vr) V 4 ) (yip + mP-Q- Q)VSlAv2P + Q, mP + Q) , (4.11) 

where y\ = (r/ + ^)/2, y 2 = (rj — £)/2 and likewise for y 1; y 2 . The factor of 1/4 on the third 
line corrects for double counting after integration over £ and £. We may approximate the 
delta-function as 

^ 4 )(y 1 p + y 1 p-Q-Q)«% 1 p+-Q+) ( 5(y 1 p--0~M (2) (Q± + 0±) , (4.12) 

which decouples the momentum fraction integrals of the upper and lower sections of the 
diagram. In addition, if we keep only terms producing a power divergence in Qj_ and use 
the gauge invariance of the Higgs vertex, Vh, we may make the replacements 

Tr fr-yw] ~ ^(2/2 - yi)^f- . (4-13) 



Tr[ 7 V^«]«4(fc-ift)^ 



P X2 



(4.14) 



Furthermore, we may write 

= EE fia2{ «- £ W • ( 4 - 15 ) 

This is the origin of the requirement that the gluons fusing to produce the Higgs have equal 



helicities, as in equation (2.3). 



Collecting everything together, we have 



2/1,2/2 



A2( - i) (4.16, 



xix 2 (iV 2 - 1) 
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Including the other diagrams, [7(b) - 7(d) and incoming gluons, we find 



dQ 2 



K ga 



l 

~2 
1 
~2 
vr 3 2 2 ( 



a,a 



da / de^a(c,^Qi)^a(e,r?;Qj 



K 9 ' a 



2 
2 



2/1,2/2 



2/1,2/2 



K 9a 



,0 



,0 



2/1,2/2 



2/1,2/2) 



xix 2 (iV 2 



(4.17) 



Here, the X denote the unregularised splitting kernels, i.e. the kernels of equations (^ 
but without the plus-prescription. The absence of the plus-prescription is due to a mis- 
match between the corrections which have the form of a self-energy and those involving 
the exchange of a parton in the s-channel. The mismatch occurs because diagrams, such 
as the one shown in figure 8(b), not involving at least one s-channel parton collinear to 
each hadron, are suppressed by the centre-of-mass energy. For now, we shall simply replace 
them with the regularised kernels but this point requires a proper treatment of soft gluon 
effects, which we shall cover in the next section. 





(a) 



(b) 



Figure 8: Two diagrams contributing to CEP. The left-hand diagram makes a leading contribution, 
whereas the right-hand diagram is suppressed by the centre-of-mass energy. 

The pdfs satisfy an evolution equation, analagous to the DGLAP equation in non- 
skewed kinematics [18]: 







din fip 
where 

and 



a s 
27 



H 



K 



(0) 



(4.18) 



(4.19) 



(4.20) 
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We may then use this equation and the symmetry relation, H g (—£, 77; Q 2 ± ) = H g (£, rj; Q\) [18], 
to write 



-4cep 



dQ 2 ± 8 
Q\ d\n{Q\) 



11 >i ( y. y;Qi 



8 



dln(Qj 



Hg[ Y , Y ,Q ± 



7T 3 2 4 (-i) - 

X (1V2 \\ V H- 

XiX2{N z — 1) 



(4.21) 



4.1.2 Corrections that generate the Sudakov factor 




Although we have now dealt with all emissions that factorise into the pdfs, there 
are still large logarithms contained in Vh- Again, we deal first with the hard collinear 
logarithms, deferring the treatment of the soft region to section 4.2. 

The diagrams contributing a collinear logarithm with respect to k 2 are shown in fig- 



ure ||, with equivalent diagrams for q collinear to k\ not shown. Dealing with diagram 3(a) 
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first, we obtain 



V, 



ab 



d 4 q H) H) 
(27r) 4 qf + ieql + ie 

The three-gluon vertex is 



f cbd V^{-q 2 M,-qi)V'^ 



acd 



(4.22) 



^^(-92,^2, -qi) 



g(9au(-q2 - k 2 )p + g v p(k 2 + qi) a + gp a {q 2 q\)v) 
- g(-(l + zi)g au k 2l 3 + (1 + z 2 )g/3 u k 2a ) 



(J 



l + Z l 1 + Z2 

-gauqi/3 + z g/3uq2a 



- Zl 1 — -22 

where we introduced the following Sudakov decomposition of the qi 



(1 - Z;)&2 + + &J- 



with 



^2 



^2 



^2 



. 



(4.23) 



(4.24) 



(4.25) 



2k 2 ■ v 

In equation (4.23) we first used the fact that we may set qf ~ (1 — Zi)k% in the numerator, 
up to terms that do not generate a logarithm. Following this, we used the transversality of 
the incoming gluons to set k% — > 0. We also multiplied and divided by 1 — Z{. This last step 
is perfectly valid for hard collinear emissions, however, it will require modification when 
we come to consider soft effects. 

Now observe that qipVlf a,acd and q 2a Vlf a ' acd describe a set of diagrams with an ex- 
ternal gluon attached whose polarisation vector is in the direction of its momentum, i.e. 
it is longitudinally polarised. The QCD Ward identity may then be applied to factor this 
gluon from V', onto an eikonal line (see for example [37]). This is shown in figure 9(a)| and 
gives the following expression 



ab 



/ Jf^lp(a) 

where 



d 4 g H) H) 

(2vr) 4 k£ + ieq 2 + ie 



f M {-gg«v{i + z)k 2 p) 



(-q-v) 



igv P { _ if dec )Vi 



ae 

Hfia 



q» = (l- z)k% + + q» 



Now, doing the j3 integral by contour integration, we find 



\rao II 1 



a s 

4-7T 



C A / dz 



(l-z)zQ 2 . 



q i- 



Vi 



Hjiv i 



(4.26) 
(4.27) 

(4.28) 



with if = q\ — (1 — z)zk\ and we made the replacement k\ ~ — Q\- Note that, since we 
are only concerned with the logarithmic terms, we may replace the lower limit on the l\ 
integral as (1 — z)zQ\ — > (1 — z)Q\. 

Before discussing the contribution from figure |9(b) , we note that these logarithmic 
corrections are going to generate the Sudakov factors, however, observe that the transverse 
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momentum integral extends down to (1 — z)Q 2 j _. This is in contrast to the integral in the 
Sudakov factor, which is cut off at Q\, though note that the difference is only relevant for 
the (1 — z)~ l term coming from diagram 3(a) . It appears then that, if we cut off the z- 
integral at 1 — z ~ \q±\/mH ( see section 4.2), this piece will generate a Sudakov factor with 
twice the double logarithmic contribution of equation (|2.7|). This is not the case however. 



Working in the transverse momentum ordered approximation (see equation fl4.4| )), we 
have missed a contribution from diagrams like those shown in figure 10(a) and |10(b) 



Clearly the diagrams in figure ^ cannot account for the screening of long wavelength emis- 
sions which results from the colour neutrality of the i-channel exchange (since they are 
oblivious to the role of the screening gluon). Fortunately, the contribution from diagrams 
like those of figure 10, in the region q\ < Q\, with the (1 — z) _1 piece generating a log- 



arithm, is included in the BFKL corrections. As we shall now see, the BFKL equation 
guarantees that there is in fact no large logarithm generated by the transverse momentum 
integral in this region and fixes the lower limit on the if integral in equation ( f4.28| ) at Q 2 ^ 
rather than (1 — z)Q 2 L . 





Figure 10: An example of diagrams generating a logarithm not included in the transverse mo- 
mentum ordered approximation. 



Consider the one- loop corrections to C qq in the BFKL region, displayed in figure 11 
where we calculate the amplitude from cuts as described in section [3|. In this region, the 
full set of virtual corrections to either side of the cut are summed up in the reggeised gluon 
propagator (denoted by a slash). These corrections amount to the replacement, in the 
lowest order graphs, of 



Q 2 (x lPl - Q) 2 (x 2 p 2 + Q) 2 Q 2 (xipi - Q) 2 (x 2 p2 + Q) 2 
where the gluon Regge trajectory, €g(Q 2 ), is given by 

Ca^s f ,d-2, Q± 



Q] 



(27T) d ~ 2 

2C A a s 
(2ir) d - 2 

C A a s (4tt) £ 
2tt r(l - e) 



d d ~ 2 k ± 



k 2 ± (k A 



d d ~ 2 k 



Ql) 2 

Q] 



(fej 



Q ± ) 2 [kl + (k ± -Q ± ) 2 } 



(4.29) 



(4.30) 
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/A. (To 



Figure 11: Next-to- leading order corrections to quark-quark central exclusive production in the 
BFKL formalism. Slashed gluon propagators and vertices labelled T indicate reggeised gluons and 
Lipatov vertices respectively, see [25]. Not shown are diagrams with the Higgs to the left of the cut. 



Note that there is no double logarithm. Moreover, the region k 2 ^ < Q\ generates only single 
logarithms in s/Q 2 ^, i.e. it does not lead to any logarithms in m^/Q 2 ^. The complementary 
region, k 2 ^ > Q 2 ^ does not generate a transverse momentum logarithm and as such it is 
never able to generate a logarithm in the Higgs mass (although it does contribute to a 
logarithm in s, which can be absorbed into the pdfs). Hence we are justified in fixing the 
infra-red cutoff in equation ( [4.28 ) equal to Q\. It will be instructive to observe how the 
elimination of the would-be double logs and the corresponding gluon reggeisation comes 
about in section ||, where we present an explicit calculation of the next-to-leading order 
corrections to C qq . 



We will close this discussion on BFKL by showing that the would-be infra-red diver- 
gence arising from ec(Q 2 ) is cancelled by emissions across the cut, which can be written 
in terms of the Lipatov vertex, VL, (see [25] equation (3.11)). After re-summation, the 
BFKL corrections can be absorbed into a function, f(oo,k±,Q±), which is related to the 
Green function with four off-shell gluons (see [25] equation (4.8)) 3 . The amplitude may 
then be computed after convolution with some non-perturbative function containing the 
long-distance physics associated with the external hadrons, as depicted schematically in 
figure Il2[ f(ui,k±,Q±) obeys an integral equation, the BFKL equation, shown diagram- 



u is a Mellin transform variable, conjugate to the centre-of-mass energy. 
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Figure 12: Schematic form of the central exclusive production amplitude in the BFKL formalism. 
The functions / are related to the four-gluon Green function and <& are proton impact factors, 
see [25]. 



matically in figure 13 



i.e. 



J (k' ± -Q ± y 

+ 2e G {-Ql)f{u,Q ± ,k ± ) 

s {2) (Q±-k ± ) 



+ 



7T 2 J {k> - Q_ 



f(u,k' ± ,k_ 



Q±f(",Q±,k±) 

Wl + (k'-Q ± y]J ■ 



(4.31) 



Now observe that, for k f <C the first and second terms in parenthesis (the first term 
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corresponds to emissions across the cut) cancel one another. 

The above is essentially the same argument, though presented in slightly different 
terms, that the Durham group use to set the lower limit on the Sudakov factor in equa- 
tion ( |2.10D . We will verify this argument explicitly in section ||, where we detail the results 
of a full next-to-leading order calculation of the virtual corrections which contribute to the 
Sudakov factor. 

Returning to the calculation of diagram 9(a)| , we find 

In addition to replacing the lower limit, we have now made explicit the upper limit of the 
l\ integration. The variable I2 is the analogue of l±, but for the next emission contained in 
V H . 

Next, the one-loop gluon propagator corrections (see for example [42]) give, for fig- 
ure 9(b) , 



vab , a s (10 2T F n f \ f^dlf b 

v H^m ~ [t 2 Ca - -3- ) j_ kl -q v w ' (433) 

where fiR is the renormalisation scale. By setting [j? R = Z|; the two contributions in figure ^ 
may be rewritten in terms of the DGLAP splitting functions, with the additional effect 
that the strong coupling now runs with if: 

f l z d/f a s (lf) f 1 / 

Vn\n collinear emissions = 2 / — ^ - / dz I zPgg(z) + TlfPqg^z) 

JQ\ h 47r Jo v 

X Vh \n— 1 collinear emissions (4.34) 

and we also included the contribution from emissions collinear to k\. The tilde again 
denotes the splitting function without the plus-prescription. 

Equation ( 4.34 ) is still incorrect as it stands, since the integral diverges in the soft limit, 
z — > 1. This is an artifact of our hard collinear approximation: for fixed, finite, transverse 
momentum, a soft gluon's energy cannot vanish. We shall discuss how this problem is 
rectified in the next section. For now, we ignore the divergence and iterate ( [4.34 ) until we 
are left with the tree level vertex: 

v* = n ii i - lT f^Bj ldz + n ^ 

dl? aM) 




f^ 1 J Q &z(zP gg {z)+n f P qg {z) 




(4.35) 

where = Z? +1 = m 2 H and M. is as defined in equation fl2.3|) . We may then use the 
following identity for ordered integrals 

/ Ml- <=\l[ (4-36) 

JQI JQl nl fA J Ql 
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which allows us to rewrite equation ( 4.35Q in terms of exponentials: 

, dl 2 a s (l 2 ' 

Vh = exp 



x exp 



Q\ 1 ^ JO 



HdPasiP) 

I 2 4tt 



dz (zPgg(z) + nfP qg (z 



dz (zP gg (z) + n f P qg (z) 



(4.37) 



x M(gg i7) . 
Thus we begin to see the emergence of the Sudakov factors. 
4.2 Soft emissions 

So far we have deferred the treatment of soft gluon effects: In section |4.1.l| we noted that 
the final rung before the hard scatter involved an unregulated splitting function, e.g. see 
equation fl4.17| ) and in section 4.1.2 we noted a potential divergence in the z — > 1 limit 



of the Sudakov factor, e.g. see equation ( 4.35 ). We shall now turn our attention to these 
matters and consider the soft gluon limit more carefully than hitherto. 

4.2.1 The Sudakov factor 




Figure 14: Form of the softest gluon attachment to the Higgs vertex. 

We shall begin our discussion of soft effects by dealing with the z — >■ 1 divergence in the 
Higgs vertex corrections. The softest gluon must attach to both gluons fusing to produce 
the Higgs, as shown in figure O. This contribution has the form 



V H \n 



soft emissions 



iC A g 2 m 2 H i v , 
(2vr) 4 



where 



d z q ± J da yd/3 

V H \ n—1 soft emissions 
^ l 

[q 2 + ie] [(q - h) 2 + is] [(q + k 2 ) 2 + is] 
q n = ak^ + + q^ . 



(4.38) 



(4.39) 
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This diagram contains divergences both when q is collinear to k\ and when it is collinear to 
&2- In order to bring the result into the form of equation Q4.34 ) we must separate out these 



two regions. An effective way to accomplish this is to multiply the integrand of ( f4.38| ) by 

1 = pv (^) +pv (^) < 440 > 

where PV denotes the Cauchy principal value. Now the term proportional to a./ (a. + f3) 
possesses only a collinear divergence with respect to k\, whereas the other piece has only 
a collinear divergence when q oc hi. We may then perform the /3(a) integrals in the 
first(second) term by contour integration. Only the piece coming from the [q 2 + ie] pole is 
relevant to the (1 — z) -1 divergence. Keeping just this piece, we obtain 



' H | n soft emissions 2, ' ™^ 



■^-3-+ Cm — ^ — s — i 



n— 1 soft emissions ■ 

(4.41) 

Changing variables as a = 1 — z, \/3\ = 1 — z and noting that the momentum fraction 
integrals are effectively cutoff at \q±\/rriH, this becomes 



Vh In soft emissions 2— Ca / o - / ~. Vff n— 1 soft emissions (4.42) 

27r J Q± Jo 1 - z 

and so equation ( [4.37 ) becomes 

V H = T(Q ± ,m H )M(gg -»• H) , (4.43) 



with 



/ r m H do 2 n (a 2 \ r l -\l±-\/ m H \ 
T(Q ± ,m H )=ex V [- -^^l dz[zP gg (z) + n f P qg (z)]\ . (4.44) 

Note that this result differs from equations d2.7|) - (|2.9| ) . In sections ^ and ^ we shall provide 



further evidence that equation ( 4.44 ) is indeed the correct form of the Sudakov factor. 



4.2.2 The Sudakov derivative 

We now turn our attention to a proper treatment of the (unregularised) splitting kernels 
entering due to the final two s-channel emissions. The only issue is with the K" kernel, 
which diverges in the soft limit. In order to understand how to treat this region correctly we 
focus on the final emission collinear to the upper hadron, assuming all previous emissions 
have been collected into the pdf. 

As discussed in section ^, we may compute the amplitude by taking the cuts shown 
in figures |3(a) and |3(e)| only. Furthermore, if we wish to study the soft limit of the final 



emission, only the diagrams shown in figure 15 and the analogous diagrams with the Higgs 



emitted on the other side of the cut will contribute. In addition, we need only consider 
the sub-set of these diagrams in which Q attaches to on-shell particles. The result is then 
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given by the soft insertion rules (see for example [43]). After summing the two diagrams 
in figure |l~5|, for example, the amplitude to the left of the cut, \Li) fl,a , may be written 

where Iq = yip, l n +\ = yip and the other ij are the momenta of the particles crossing 
the cut. The amplitude without the soft gluon is represented by \Lo). Since Q is the 
last emission collinear to the upper proton, we may take li oc p for i / 0. Including the 
amplitude to the right of the cut, \Ri), using the same soft insertion formula and including 
an integration over the intermediate phase-space, we obtain 

n+1 

£ ((Ro\ {T^Tt |L ) + <i2o| (TtjTH \L )) , (4.46) 



n+1 

X 
i=l 



where d(PS n ) is the phase-space of the cut diagram without the soft gluon. Then, using 
colour conservation: 

n+1 

E T £ = - T £> ( 4 - 47 ) 

i=l 

{T-)%1 = C A , (4.48) 

we find 

Now note that (|4.49| ) possesses divergences when Q becomes collinear to either p or p. As 
in the case of the Sudakov factor, we may separate out these regions by multiplying the 
integrand by 1 = (a + /3)/(a + /3), where this time 

Q» = ap» + ffi + . (4.50) 




Figure 15: Diagrams for the last emission collinear to the upper hadron which generate the soft 
part of the splitting kernel. 
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Since we are interested in the divergences with respect to the upper hadron, we keep only 
the a I (a + f3) piece. The only effect of the phase-space on Q is to introduce the constraint 
®(yiP + ~Q + )- Thus, flOp becomes 



A 



2g 2 C A f d 2 Q 
(2^) 3 J Ql 



da- 



a + 



91 

Acts 



d(PS n ) (R \L C 



(4.51) 



Comparing this to equation ( |4.5| ), by taking the limit Xi — > yi, we see that the correct form 
of the unregularised splitting kernel is 



Xl Xl ^(xijZi - 2/1) 



+ 



(sci— j/i)4s y 

x 2 x 2 •&i 1 (x 2 ,X2 - 2/2) 



2/2 (#2 - 2/2 + 
+ 2 XlX2 + /^ 2 ^"11 (ai, -aJ2,ai - 2/1) 

- 2 "- — , ToX ^!, -^2) 



, xix 2 sij/i + x 2 y 2 nQ p 
2/12/2 (xi + x 2 ) 2 



(4.52) 



To the soft divergent pieces we may now add and subtract a term proportional to a delta- 
function: 



x x $ii(x, x — y) 
y [x - y + 



Q 1 



(x-y)4s 



X 

y 



x $ii{x, x — y) 



x - y + 



Q 



(x-y)4s 



-5(x — y) dx 



, x'$\ x (x\x' -y) 



y + 



Qj 



+ S(x — y) dx 



(x'-y)4s 

, x'$° n (x',x' -y) 



1 , gl 

» ~ (x'-y)4s 



(4.53) 



In the first term, contained in square brackets, we may take the Q± — > limit since the 
soft region cancels, leaving us with something regularised by the plus-prescription. The 
integral in the second term may be done explicitly, giving a logarithm plus non-logarithmic 
terms which we neglect. Substituting this identity into equation ( 4.52j ) it becomes 

'Ayls + Ql 



K"(x 1 ,x 2 \y 1 ,y 2 ) = K"(x 1 ,x 2 \y 1 ,y 2 ) - —C A S{xi - y\) In 

47T 



Q 2 



a 



—C A S(x 2 -y 2 ) In 

4-7T 



42/jg + Q] 

Ql 



Applying this to the splitting kernels appearing in equation ( [4.17 ) gives, e.g. 

//1-//2 I - <>'"—( ' i<> K I In 



K° a 



2/1,2/2 



K<> a 



m H 

Ql 



(4.54) 



(4.55) 
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where we replaced x\s — > m 2 H , which is correct to logarithmic accuracy. This logarithm in 
the Higgs mass may then be written in terms of the derivative of the Sudakov factor: 



K^ a ['^, o 



) A'"" (y,0 
- 5 9a 5 U 



2/1)2/2 



dy/T(Q ± ,m H ) 



2J VT(Q±,m H ) dln(Ql) 
All together then, including the Sudakov factor, equation ( }4.2l| ) becomes 



(4.56) 



-4.CEP 



dQ 







Qi dHQl) 







dHQl) 

vr 3 2 4 (-i) 



X2 X2 

T'T 



-M(gg^H) . 



(4.57) 



xix 2 {N 2 - 1) 

Finally, using equation fl2.6| ) and assuming that R g depends only weakly on Qj_, we find 

fdQ 2 

Acep ~ J -^fg(xi,0,Q±,m 2 H )f g (x 2 ,0,Q\,m 2 H ) 

7T 3 2 4 (-i) 



x lX2 {N 2 - 1) 



M(gg^H) . 



(4.58) 



5. An explicit next-to-leading order calculation 



Having discussed the result at all orders, we now turn to a description of our next-to- 
leading order calculation of the amplitude for two quarks of different flavour to scatter into 
two quarks and a Higgs. This calculation will serve as an explicit check of the all orders 
result presented in the previous section and also offers the possibility to extend the result 
to next-to-leading order accuracy 4 . 

As stated in section |3|, we may calculate the amplitude from the cuts 3(a), 3(e). We 
limit ourselves to a calculation of the virtual corrections to one side of a cut; it is the 
Sudakov factor we are interested in probing here and we expect only this set of diagrams 
to contribute to it. After presenting our results we shall comment on why the diagrams 
with an additional gluon crossing the cut cannot contribute to the Sudakov factor. 

The set of diagrams we must calculate when the Higgs is to the right of the cut are 



shown in figure 16, with a similar set for the Higgs to the left of the cut not shown. All 
other diagrams may be obtained by exchanging x\ and x%- We perform the loop integrals 
using the techniques described in [46] , which we have implemented using Mathematica [47] 
and FORM [48]. We also use the Mathematica package FeynCalc [49] to simplify the 
numerator algebra. All of our calculations are performed in Feynman gauge and using the 
large top mass effective theory described in appendix [A]. 



4 The Durham group do include a K-factor in their calculation of central exclusive Higgs production [9], 
taken from the calculation of inclusive Higgs production at next-to-leading order [44, 45] . They do not 
however explicitly evaluate the next-to-leading order contribution. 
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In addition, we keep only terms not suppressed by additional powers of Xi or Q 2 ± , 
relative to the lowest order case. The cut sets 

Q ± - ^ (5-1) 
and we may also make the approximation 

UaM ~ u ai {Pi) , (5.2) 
for the final-state quark spinor wavefunctions. This allows us to make the replacement 

UaMua^p'i) , (5.3) 

where we have dropped a term proportional to 7 5 , which is not relevant since the amplitude 
is CP- invariant. We shall also present here only those terms either enhanced by a logarithm 
or divergent as e — > 0. 

Working in dimensional regularisation, with d = 4 — 2e space-time dimensions, the 
results for each cut diagram are 
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(5-4) 
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CnL0||16(s1 = CNLo||l6(d) 



(5.20) 



C 



NLOll6(t) 
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NLO|l6(c) 



(5.21) 



Cnlo||i6(^ 



CNLo||l6(f) 

where the usual MS factor is given by 

M = exp [e(-j E + ln(4vr))] 



(5.22) 



(5.23) 



and je is the Euler-Mascheroni constant. The full set of counter terms for the diagrams 
with the Higgs to the right of the cut give 



Cnlo 



counter terms 



dQj 
<51 



IT € 



where /3q is the first component of the QCD beta-function and is given by 

11C A -AT F n f 



A) 



12 



(5.24) 



(5.25) 



Collecting these results together and including those related by interchanging x\ -H- X2 
and those generated by the diagrams with the Higgs to the left of the cut, we obtain 



-4nlo = ^o(m) J 



dQj 



a s M ( n 2 
vr e 2 \Q\ 
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(5.26) 



The 1/e poles here are due to collinear and soft divergences and may be written in terms 
of the quark-quark splitting function using the following identity 



C F 



(Qir e , 3(Qir 



+ 



Q 



2 \l+e 



1-9x/IQj 



dz P qq (z) , 



(5.27) 



which holds up to terms not involving either 1/e poles or logarithms of Q±. Each of 
these factors, proportional to P qq , is due to a parton becoming collinear with one of the 
on-mass-shell, external quark lines. The (final state) contributions associated with the cut 



- 28 - 



quark lines will cancel with an equal and opposite term coming from final-state divergences 
associated with a gluon emitted across the cut. If we take fi 2 = Q\, then the remaining 
factors of P qq are simply the virtual contribution to the expansion of the pdf, H q (£, 77; Q\)- 
Finally, we see the factor associated with the reggeisation of the gluon, written in terms of 
the gluon Regge trajectory (see equations J4.29 ) and ( |4,30[ )). This has precisely the form 
predicted in equation ( f4.29| ) and is also expected to factorise into the unintegrated gluon 
pdf. 

What remains then, after accounting for all of these pieces, must be the 0(a s ) expan- 
sion of the Sudakov factor: 
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It is interesting to see which diagrams this double logarithm derives from. If we simply 



took the vertex correction of figure 16(a) , we would obtain twice the double logarithmic 
contribution of the full answer, which is consistent with our discussion of the corrections 



to Vh in section 4.1. This additional double logarithm is cancelled when we include the 
diagrams 16(g), 16(h) and 16(m)- 16(p)| (and those related by x\ o £2). It is also the sum 
of these diagrams which generates the reggeisation of the gluons to the right of the cut. 

In order to compare equation ( |5.28 ) to the Durham result, we must first exchange 
Aq(Q±) for Ao(mH)- This may be accomplished using the leading order coupling constant 
evolution 



a s (/4) = a s {nl) exp 



which implies at next-to-leading order 
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This is exactly the next-to-leading order term one would obtain by expanding out the 
Sudakov factor (with A = k^jrau) to this order in perturbation theory. 

In obtaining this result we have neglected the diagrams in which a gluon is emitted 
across the cut. However, we now argue that such terms cannot possibly contribute to the 
Sudakov factor. 

There are three regions for these diagrams which could potentially produce a loga- 
rithm in m'j]/Q 2 L . First are the BFKL corrections. However, these are summed into the 
unintegrated pdfs and so cannot contribute to the Sudakov factor. Secondly, we have hard 
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collinear emission. This was fully accounted for in section 4.1 , where we saw that the only 
large logarithms generated by s-channel emissions in this region had transverse momentum 
less than Q± and so do not form part of the Sudakov factor. Finally then, we have soft 
emission. To be included in the Sudakov factor, the additional soft gluon must have a 
larger transverse momentum than one of the quarks. However, this means it is one of the 
final two s-channel emissions making up C| s-channel ■ ^ s snowri i n section 42, these soft 
emissions generate the Sudakov derivative. 



With these considerations then, we conclude that equation ( 5.31|) gives the full con- 
tribution to the Sudakov factor at this order in perturbation theory. Moreover, we have 
confirmed our previous result for the resummed Sudakov factor, equation ( 4.44| ). 
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Figure 16: Virtual corrections contributing to the cut qq 1 — > q © H © q 1 amplitude at next-to- 
leading order, in the high energy limit. Not shown are those diagrams obtained by exchanging x\ 
and X2 and those in which the Higgs is radiated on the left of the cut. 
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(a) (b) (c) 



(d) 

Figure 17: Diagrams contributing to the process gg —> Hg. 
6. A further cross-check 

We now discuss a re-calculation of the soft part of the Sudakov factor, which exploits 
the Bloch-Nordsieck theorem [50]. It serves as a further check on our earlier results. In 
particular, we infer the virtual corrections by considering the set of real emission diagrams 



making up the gg — > gH amplitude, as shown in figure 17. The relevant cross-section is 
given by 

v = ^6 J d 6 s(£i;/40s(&;/4) f da f d/3 f J2~^ 

x (27r)5 (+) ((fci + k 2 - q) 2 - m 2 H )\M\ 2 (6.1) 

where the & denote the momentum fractions of the incoming gluons and we parametrise 
the final-state gluon momentum, q, in terms of Sudakov variables as in equation ( 4.3S| ), 
with 

s = 2fci-fc 2 - (6.2) 

The amplitude, Ai, is summed over equal incoming helicities by contraction with a polar- 
isation vector (in the k\-h,2 centre-of-mass frame), = (0, cos </>, sin <j), 0) for the incoming 
gluons, followed by an average over <f> i.e. 

M = p-e^M^ • (6.3) 
Jo 27r 

Next, we divide the phase-space of the emitted gluon into the regions a > j3 and a < f3, 
which is equivalent to dividing the phase-space about zero rapidity. We shall calculate only 
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the contribution from the a > (3 region, however it will be clear that the a < (3 region 
gives an equal result. In addition, we sum over final state polarisations using 5 



fci • k 



(6.4) 



We also average over incoming gluon colours, though we stress that changing the treatment 
of colour in the process affects only the normalisation, since each diagram has identical 
colour structure. After a little algebra, which is performed using FORM, we find 
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(6.6) 



is the lowest order gg ^ H amplitude summed over equal incoming helicities. The key 
point to note is that only the term proportional to \/q\ can generate a logarithm or a 
constant after the a integral: The remaining terms are insufficiently singular in the soft 
(a — > 0) limit. Keeping only this term then, changing variables as a = 1 — z, approximating 
the delta- function and including the region a < (3, we find 
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(6.7) 



where we replaced y§ — > niH in the upper limit of the z integral, as appropriate in the 
soft limit. The additional z-dependent piece, not proportional to P gg {z), arises due to the 
restricted sum over incoming helicities and it vanishes in the soft z — V 1 limit. 6 Thus, in 
the soft limit, we find 



C soft OC 



dq\ C A a s 



In 



?i f \q± 

which further supports the validity of equation ( 4.44j ). 



(6.8) 



Note that, in processes involving two or more external gluons, the replacement of the gluon polarisation 
sum with — <7 M " is not in general valid [51]. 

6 If we had summed over all incoming helicity configurations this term would of course have been absent. 
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Figure 18: Ratio of the cross-section for central exclusive Higgs production at the LHC evaluated 
with the scale in the Sudakov factor set to \i — m# divided by the cross-section with the scale set 
to fi — 0.62mjy, plotted as a function of the Higgs mass. The solid blue and dashed red lines were 
generated using MRST2002nlo [53] and CTEQ6m [54] parton distributions respectively. 



7. Cross-section predictions 

We can assess the impact that our modification of the Sudakov factor has on predictions 
of the central exclusive cross-section. Taking, as an example, the cross-section for central 
exclusive Higgs production at the LHC, with 14 TeV centre-of-mass energy we compute 
the cross-section, using the ExHuME Monte Carlo generator [52], placing no cuts on the 



final-state particles. The results are shown in figure 18, for two different parton distribution 
functions [53,54]. We observe a suppression of the cross-section, relative to the predictions 
of the Durham group, by a factor ~ 2 which increases with increasing Higgs mass. 



8. Conclusions 



In this paper we have studied the central exclusive production process in some detail. In 
particular, we studied the cross-section for Higgs boson production, using QCD perturba- 
tion theory. We largely confirm the calculation previously performed by the Durham group 
(e.g. see [27]), except that we disagree as to the precise form of the Sudakov factor which 
enters. Using the Sudakov factor that we propose leads to a suppression of the central 
exclusive production cross-section at the LHC by approximately a factor of two relative to 
the earlier predictions, for Higgs boson masses in the range 100-500 GeV. 

As a point of further study, it would be interesting to assess the impact on predictions 
for other processes and in particular on the central exclusive production of dijets at the 
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Tevatron, for which data exist. We do not expect to find any disagreement with the data. 
In fact, the reduced suppression at lower masses suggests that agreement with the data 
may even be slightly improved. However, one must always remember that the theoretical 
uncertainty on other parts of the calculation, for example the gap survival factor and 
unintegrated pdfs, is expected to be comparable in size to the effect induced by the change 
in the Sudakov factor that we have been focussing on (e.g. see [55]). 

We note that the fixed-order corrections we have computed form a subset of the full 
next-to-leading order corrections to central exclusive Higgs production, offering the possi- 
bility of extending the theoretical description of the process to this order. 
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A. Large top mass effective theory 

In this appendix we describe the effective theory, formed by taking the top quark mass to 
infinity [28-30] which we use to compute the next-to-leading order corrections to central 
exclusive Higgs production. This approach has been found to give good agreement with 
the full theory, provided that the Higgs mass satisfies mn < 2mt op and the transverse 
momenta of any jets produced in association with the Higgs satisfy p± < rn top [56,57]. 

We work in a theory in which the top quark has been integrated out and all other quarks 
are taken as massless. In this approach, the only coupling of the Higgs is to gluons 7 , via 
the following term in the effective Lagrangian 

C eS = -j C?0? , O? = (Go)%(G Q r» (A.l) 

where Cf is a coefficient function, zeroes indicate that these are bare quantities and it is 
understood that they are defined in the five flavour effective theory. Both C® and matrix 
elements of 0\ contain ultraviolet divergences, however their product is finite, since the 
operator in the full theory which ( |A.l| ) approximates (^mtV?^) is a conserved current. 

Using the Bogolyubov-Parasiuk R-operation [58,59], it is possible to define a finite 
version of O®, which we denote O^. This finite operator may then be related to the bare 
operator as [59] 

Of = Z Ol O l , Z 0l = 1 (A.2) 

1 - p(Na 8 )/e 

where 

AA = exp[e(- 7E + ln(4^))], (A.3) 



7 We consistently ignore electroweak couplings throughout. 
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Figure 19: Feynman rules for the large top mass effective theory. See the text for the definitions 
of V3 and V4. 



and f3{a s ) and a s = a s (fi) are the QCD beta function and the MS running coupling 
respectively. Again, both are defined in the five flavour theory. 
The effective Lagrangian now reads 



H 



-cff 



R _ ry-Xrfo 



(A.4) 



Since both matrix elements of Of and the full expression are finite, Cf is also finite and 
is given by 

1 a s (ji) ( 11 a s {\i) 



3t> 7T 



1 + 



4 7T 



+ OK 



(A.5) 



where again a s (fi) is the five flavour MS running coupling and v is the Higgs vacuum 
expectation value. 

The Feynman rules generated by equation (|A.4 ) are displayed in figure 19, and can be 
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written in terms of the standard three and four gluon vertices: 



^^i^mIm^ 1 ' ^ 3 > ^4) = — ^ (f aia2e f aia3e (9mn4,9n2^3 ~ 9fnn3 9^4^) 

fa ia , 3 e ra 4 a 2 e / _ \ 



(A.6) 
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